Detection of transverse entanglement in phase space 
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Transverse entanglement between pairs of photons can be detected through intensity correlation 
measurements in the near and far fields. We show theoretically and experimentally that at in- 
termediate zones, it is also possible to detect transverse entanglement performing only intensity 
correlation measurements. Our results are applicable to a number of physical systems. 



Detection and quantification of entanglement is es- 
sential for the development of many applications in 
the field of quantum information. Several tasks pro- 
posed to take advantage of the entanglement properties 
of quantum systems can be experimentally tested with 
photons produced from spontaneous parametric down- 
conversion(SPDC) fl. This is a versatile system, since 
SPDC photons can be prepared in entangled states of 
many different degrees of freedom, such as polarization 
0|, time-bins 0], orbital angular momentum 0], as well 
as transverse spatial variables 0, Q . The latter concerns 
correlations between the transverse components of the 
wave vectors of the signal and idler photons, which have 
been extensively studied and utilized in the last decade 
They arise due to the localization of the emis- 
sion of photon pairs and the phase matching conditions 
for the non-linear interaction between the pump, signal 
and idler fields. Even though the quantum nature of spa- 
tial correlations was already evident [Io| , the formal rela- 
tionship with entanglement has only been demonstrated 
a few years ago 0, Q . Transverse entanglement was de- 
tected through the violation of a non-separability criteria 
Ul,LL2| , based on intensity correlation measurements per- 
formed in the near and far fields. Entanglement in con- 
tinuous variables (CV) is a rich research subject, because 
several quantum information tasks can be optimized us- 
ing high dimensional Hilbert spaces 114 1 . SPDC is 
a natural option for the experimental investigation of 
transverse spatial entanglement, which can be present 
in other quantum systems (l5l . 

So far, CV transverse entanglement detection has been 
based on intensity correlation measurements performed 
in the near and in the far field [H, @]. An interesting 
entanglement "migration" effect was shown recently by 
Chan et. al (l7j . in which entanglement moves from the 
real to the imaginary part of the two-photon wave func- 
tion during propagation. In order to be able to detect 
entanglement in this case, it would be necessary to per- 
form phase-sensitive measurements [lij ]. 



In this Letter, we show theoretically and experimen- 
tally that it is always possible to detect entanglement 
by performing intensity correlation measurements, even 
outside near and far field zones. We demonstrate the 
connection between the variances of two observables and 
the variances of these same observables rotated in phase 
space. We encounter the conditions for which entangle- 
ment detection is possible with intensity measurements, 
and others for which it is impossible. This connection 
allows one to circumvent problems like the migration of 
entanglement fl7| by performing proper phase space ro- 
tations on the observables. Though wc consider the par- 
ticular case of propagation of transverse correlations of 
photon pairs, our results can be used to improve detec- 
tion of entanglement in other CV systems. 

Our approach is based on the propagation of the sig- 
nal and idler fields using the formalism of the Fractional 
Fourier Transform (FRFT), which is parameterized by 
the angle a [3] . The FRFT appears naturally in a num- 
ber of physical systems and describes rotation in phase 
space. In particular, it is possible to completely describe 
the propagation of a light field through the order a of the 
FRFT. For instance, the field at the source is given by a 
FRFT of order a = and the usual Fourier transform, 
associated to Fraunhoffcr diffraction in the far field, is 
given by an FRFT of order a — tt/2. Free propagation 
can always be described in terms of an FRFT operation 
up to a quadratic phase term, which can be considered 
essentially unity in the near and far field [24[. In Ref. 
0], as is customary, detection of entanglement was actu- 
ally performed using lenses to obtain the intensity cor- 
relations in the near field image (a = tt) and far field 
(q = 7r/2). Likewise, any FRFT of order a can be im- 
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plcmcntcd perfectly with lenses 

We consider the experimental arrangement sketched in 
Fig. [I] a), where signal and idler photons from SPDC are 
sent through FRFT systems of order a s and a>i, respec- 
tively. 

Following Ref. [13], we write the two-photon wave- 
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FIG. 1: (Color online) a) Experimental Setup. Boxes on sig- 
nal and idler paths represent the optical systems used to per- 
form FRFT of order a s and Oj. b) The optical lens system 
used to perform the FRFT [2(J. L is a lens with focal length 
"/" and z a = 2/ sin 2 (a/2). 



function at the source in position representation as 

(Pi+Ps) 2 (Pj-Ps) 2 
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where \^f) is the two photon quantum state produced by 
SPDC, p s and pi are appropriate dimensionless position 
coordinates in the source plane and <7 + and <j_ are in- 
dependent functions of experimental parameters. Once 
the state at the source is taken to be approximately gaus- 
sian, we only have to consider one spatial dimension. The 
two-photon wavefunction in wavevector representation is 



(9 s ,ftl*> 



(2) 



where q s and qi are transverse dimensionless wavevector 
components at the source plane. 

In order to detect entanglement of the state \^>), one 
must apply a separability criteria ll|, [12J. For exam- 
ple, defining the dimensionless operators p± = pi ± p s 
and q± = q { ± q s {[p Ak] = i5j,k j,k = s,i) the 
separability criteria of Duan, Giedke, Cirac and Zollcr 
(DGCZ) [ll| establishes that if one of the two inequalities 
<(Ap_) 2 }* + ((A<? + ) 2 )* > 2 or ((A / 5 + ) 2 )* + ((Aq_) 2 )* > 
2 is violated, then the state |>F) is non separable and 
therefore is entangled. 

Using ((T|) and we have 



(A0S + ) 2 )*=a 2 , 
(A(p_) 2 )*=a 2 , 
(A(g+) 2 )* = l/a 



(3) 
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and we obtain 



((A/M 2 }* + {(Aq + ) 2 h 



(4) 



The right hand side (RHS) of Eq. (|4]) can be smaller 
than 2 for small <r_ and large a+. In the case of SPDC, 
this is readily achievable, as these two parameters are 
independent and experimentally accessible. 

The variances in inequality @ refer to position and 
momentum variables of the signal and idler fields in the 



source plane, which are related to the intensity distri- 
butions in the near and far field. It is well known that 
the propagation of a light field characterized by a FRFT 
is equivalent to a rotation of the transverse variables in 
phase space given that these variables are prop- 

erly adimcnsionalized 



2l| . The dimensionless operators 



transform as 
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sin ctj qj 
- cos otj qj , 



(5) 



where j = s,i. Therefore, it is possible to write the 
DGCZ inequality for rotated transverse variables of the 
fields, p'_ = p as — p ai and q' + = q as + q ai , in terms of the 
variables p_ and q + at the source j22[: 



((Ap>_) 2 )* + ((Aq> + ) 2 )* = 
l+cosC^+a.) [ ((A p_) 2 )^+((Ag + ) 2 ) g >] 



(6) 
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[({p + ,q + })-2(p + )(q + )U 
+ ^^ [({p-,q-}}-2(p_)(q.}] 9 . 

Eq. ([5]) shows that whenever a% + a s (mod 2tt) = , the 
sum of variances for the rotated variables coincide with 
the sum of variances for the variables at the source. This 
shows that, for any propagation of the signal field, char- 
acterized by a s , it is possible to find a propagation of the 
idler field ctj, so that intensity correlation measurement 
will violate the DGCZ inequality, in or out of the near 
and far field. We also note that Eq. ([6]) does not depend 
on the state (Eq[T]) and is applicable to any bipartite 
continuous variable systems. 

For states of the form ([T]) the last two lines of the RHS 
of Eq. @ are zero. Then considering an entangled state 
satisfying ((Ap_) 2 )* + ((Aq + ) 2 )* = ct 2 + l/a\ < 2, a 
necessary condition to detect entanglement is 



Si + S 2 - 4 
cos(ai +a s )> — - S2 > 0, 



(7) 



where we define Si = a\ + 1/er 2 and S2 = cr 2 + 1/er 2 . 
We note that for cos(a.; + a s ) = 0, intensity correlation 
measurements never evidence entanglement, regardless of 
the state. 

We have experimentally tested these conditions, using 
pairs of twin photons generated by SPDC in a 5mm long 
lithium iodate crystal (LHO3) with a c.w. diode laser os- 
cillating at 405nm, as shown in FIG.Q]a). Optical FRFT 
systems, such as the one shown in FIG. [T]b) were used 
in each of the down-converted fields. This system, with 
z a = 2/ sin 2 (a/2), is able to implement a FRFT in the 
range < a < n. For a > n we use a series of FRFT 
systems, respecting the additivity condition of maintain- 
ing /' = / sin a [l9( the same for each. To describe all 
FRFTs as rotations in the same phase space, we use di- 
mensionless coordinates p = \fkJY 1 P an( l 9 = \J f'/k Q- 
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In our experimental setup, /' = 25/\/2cm. Signal and 
idler photons were detected with single photon counting 
modules and lOnm FWHM bandwidth interference fil- 
ters centered at 810nm. Horizontal slits (100pm) were 
mounted on translation stages and scanned vertically in 
steps of 50/xm to register the detection position. In all 
measurements, the "+" ("-") correlations were measured 
in all cases by scanning the detectors with equal steps in 
the same (opposite) directions. 

First, we measured the p_ and g+ distributions at the 
source, using imaging (a a = a, = ir) and Fourier trans- 
form (a s = on = 7r/2) lens configurations The di- 
mensionless variances were A 2 (p_) = 0.93 ± 0.01 and 
A 2 {q+) = 0.073 ± 0.004. Applying the DGCZ inequality 
we obtain 



a) 
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A(p_) 2 + A( g+ ) 2 = 1.00 ± 0.01 < 2, 



(8) 



indicating that the state is entangled. 

Next, we measured the intensity correlations for the 
signal and idler fields at intermediate zones. We chose 
FRFT orders a s — on = 37r/4, so that cos(a s + a,-) = 0, 
which does not satisfy the condition of Eq.Q. The co- 
incidence counts C(p\„ — pf^) and C(p l s 2L + p%- K ) are 

4 4 4 4 

plotted in FIG. |2J-a) and (J2]-b), respectively We obtain 
A 2 (pL - PL) = 13.6443 > 2 and A 2 {p^ + p%„) = 

4 4 4 4 

39.1473 > 2, which clearly indicates that these intensity 
correlations cannot be used to violate the DGCZ inequal- 
ity. We also tested an intermediate zone configuration 
following the condition given by Eq. ([7]). We used three 
additive FRFT lens systems to perform a a s = ^f- order 
FRFT on the signal field, while maintaining the a 
order FRFT on the idler field, so that on + a s 
Coincidence counts C(p l 3 JL — pf^) are plotted in FIG. [2} 

4 4 

c), and the dimcnsionless variance is A 2 (p l ^ JL — plU) = 

4 4 

0.038 ±0.005. Coincidence counts C(gL ±gL), plotted 

4 4 

in FIG.J2J1), were measured performing an inverse Fourier 
transform of the signal and idler fields at the planes of 
FRFT of order and ^, corresponding to FRFT of or- 
ders ^ and respectively. The dimcnsionless variance 
is A 2 (q| 2L + q% L ) = 0.069 ±0.003. With our experimental 

4 4 

data, we are now able to verify entanglement at interme- 
diate zones: 



37T 

4 

2tt. 



A 2 (pL-p% L )+A 2 (q\ 



0.107±0.006 < 2. (9) 



The experimental values obtained in Eqs. (JSj) and (|9|) 
are not equal as expected from Eq. ([6]) . This discrepancy 
can be explained by the experimental imperfections. It 
is difficult to characterize every source of experimental 
error and their precise effect on the measurement results. 
However, we notice that these imperfections contribute 
by broadening the coincidence distributions. In this re- 
spect, our measurement results are upper limits to the 
actual variances. To evaluate the expected variances, we 
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FIG. 2: (Color online) Measured coincidence counts and 
gaussian curve fits, a) C(pV„ — P%k) b) C(p l 3jL + p%v_) c) 

C(pL-p^) d) C(qL+q£) 



characterized the initial state (p} at the source by mea- 
suring the width w of the intensity distribution of the 
pump beam. The dimcnsionless variance at the source 
IS <7l = (Aiu 2 )k/f = 47 ±2. The dimensionless vari- 
ance er 2 is given by cr 2 = (k/f')0A55D/K = 0.006, 
where D is the length of the nonlinear crystal and K 
is the pump beam wavenumber (l7| . With these val- 
ues, we predict a violation of the DGCZ inequality: 
((Ap-) 2 )^ + ((Aq + ) 2 ). t/l = 0.027 ± 0.001 < 2, which 
is smaller than both experimental values ([5| and ©. 
Therefore it is clear that without the experimental im- 
perfections, we should have observed even stronger vi- 
olations for both cases. A possible imperfection is the 
error in lens positioning. The effect of this type of im- 
perfection can be estimated by calculating the propa- 
gation through the different lens systems in each field 
[19( and including a 1% error in z for all lens systems. 
Taking the worst case scenario, we obtain the following 
predictions for each variance: A 2 (p l s 2L — pij-)^ = 0.09 

and A 2 ( 9 L + qU)th - 0.04, A 2 (p_) th = 0.84 and 

4 4 

A 2 (q + )th = 0.03. These variances are much closer to 
the experimental values. Thus, considering small exper- 
imental imperfections, the theoretical prediction agrees 
with both Eqs. ([8|) and ((9]), as expected from Eq. ([6|). 

Let us now discuss the application of these results to a 
situation similar to that of Ref. (l7l ] , in which it is shown 
that the transverse intensity correlations decrease as the 
field propagates, and then are recovered again in the far- 
field. For a certain propagation distance, the DGCZ or 
similar inequality will be satisfied, because the real part 
of the wavefunction becomes separable and the entangle- 
ment is present only in the imaginary part. An analysis 
similar to that of Ref. 17|] in terms of FRFT's yields a sep- 
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Pump 



describes rotation in phase space, our results are appli- 
cable to a number of physical systems. 
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FIG. 3: (Color online) Both signal and idler fields propagate 
through lens systems that implement a FRFT of order a sep = 
tan _1 (<7+<7_). The dashed detectors will observe no intensity 
correlation. Additional FRFT systems are used so that the 
total propagation for signal and idler are characterized by 



FRFT's of orders et s 



ctsep + a S 2 and an 



+ di2, 



respectively. Intensity correlation is recovered for any a S 2 
and «i2 such that a; + a s (mod27r) = 0. 



arability condition for the real part of the wavefunction 
([1]) given by a sep = tan _1 (er+cr_), where a s = on = a sep 
is the order of the FRFT implemented on both fields. 
Substituting this condition in Eq. Q shows that the 
DGCZ inequality is not violated. 

To successfully detect entanglement in this case, it is 
necessary to adopt a scheme such as the one shown in 
Fig. [3l Signal and idler fields propagate through opti- 
cal systems characterized by FRFT's of order a sep , so 
that the separability condition for the real part of the 
wavefunction is reached. At this point, intensity correla- 
tions alone will fail to register entanglement. To detect 
entanglement, additional optical systems are used to im- 
plement a second FRFT in the signal and idler beams, 
with orders a s2 and a i2 respectively. The optical systems 
can be designed so that the additivity of two consecu- 
tive FRFTs is preserved The signal beam has had 
a total propagation characterized by a FRFT of order 
a s = Q-sep + oisi and the idler ai = a sep + a>i2 ■ In this 
case, applying the condition ai + a s (mod27r) = 0, Eq. 
(|6]) shows that the DGCZ inequality will be violated in 
the same way as it would be for the field in the source. 

In conclusion, we have demonstrated theoretically and 
experimentally, that it is possible to detect transverse 
entanglement performing intensity correlation measure- 
ments, not only in the near and far fields, but also at 
intermediate propagation planes. This is achieved using 
optical systems that implement Fractional Fourier Trans- 
forms (FRFTs) according to the condition a s +cti = 2nir, 
where a s and ai are the orders of the transforms on the 
signal and idler fields. We also show that entanglement 
is never registered when a,; + a s (mod7r) = n/2. These 
results demonstrate that, given a signal field propaga- 
tion characterized by a s , one can always find an FRFT 
ai which can be used to detect entanglement with inten- 
sity correlations alone. Though our experiment was con- 
ducted using spatial entanglement of photons, our results 
are directly applicable to spatial entanglement in other 
systems 15|, [iff . Since the Fractional Fourier Transform 
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